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I.  INTRODUCTION 


Intensity  Modulated  Radiation  Therapy  (IMRT)  is  the  most  technologically 
advanced  treatment  method  available  in  external  beam  radiation  therapy.  In  recent  years 
the  advantages  of  IMRT  have  been  well  established  and,  not  surprisingly,  there  has  been 
increased  reliance  of  IMRT  accelerator-based  radiation  therapy  versus  3D  conformal 
treatment’. 

With  IMRT,  the  radiation  beam  is  effectively  broken  up  into  thousands  of  pencil- 
thin  radiation  beams.  This  approach  takes  advantage  of  steep  dose  gradients  allowing  for 
a  more  conformal  treatment  plan,  which  translates  to  improved  sparing  of  normal  tissue. 

A  common  technique  of  creating  IMRT  beams  is  to  employ  a  segmentation  algorithm, 
which  then  converts  an  ideal  fluence  map  into  deliverable  beam  segments.  These  pencil 
thin  beams  are  formed  by  multileaf  collimators  (MLC’s),  which  narrow  down  a  larger 
beam  into  a  series  of  small  rectangular  fields.  These  small  rectangular  fields  are  typically 
1.0  cm^  or  larger.  While  this  technology  represents  a  vast  improvement  over  typical 
external  beam  therapy,  it  also  may  present  a  host  of  potential  problems  when  it  comes  to 
treatment  planning. 

During  treatment  using  these  MLC  beams,  the  separation  distance  between  the 
MLC’s  may  be  as  little  as  5  mm  for  delivery  of  IMRT.  One  of  the  challenges  of  this  type 
of  delivery  is  whether  or  not  the  treatment  planning  software,  such  as  AD  AC  Pinnacle, 
can  accurately  model  these  small  fields  inherent  to  IMRT  plans  and  provide  accurate 
dosimetry  data.  Small  field  dosimetry  in  itself  presents  a  challenge.  Issues  of  penumbra, 
scatter,  and  volume  effect  can  all  adversely  affect  data.  Add  to  this  the  anomalies 
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associated  with  MLC’s  i.e.,  interleaf  leakage,  intraleaf  transmission,  leaf-end 
transmission,  and  one  begins  to  realize  how  daunting  the  dosimetric  challenges  can  be. 

Penumbra  is  the  region  at  the  edge  of  the  beam  over  which  the  dose  rate  changes 
rapidly  as  a  function  of  distance  from  the  beam  axis.  In  one  study,  Bayouth  et.  al. 
determined  the  penumbra  for  Icm^  fields  to  be  0.36cm  +/-  0.03cm  for  99%  of  the 
measurements.  Proper  penumbra  calculation  is  critical  for  IMRT  treatment  planning 
systems.  The  position  of  the  MLC  leaves  depends  on  the  penumbra  calculated  by  the 
planning  software.  Inaccurate  calculation  of  penumbra  may  result  in  cold  or  hot  spots 
between  two  adjacent  fields. 

With  regards  to  the  detector  size  (volume  effect)  in  the  dosimetry  of  small  fields, 
Laub  et.al.^  determined  that  while  using  a  0.6  cm^  Farmer  chamber,  deviations  between 
the  absolute  point  dose  values  calculated  in  ADAC  Pinnacle  and  dose  values  measured 
with  the  chamber  were  as  high  as  6%.  With  a  pinpoint  chamber,  (e.g.,  PTW-Freiburg 
diode  type),  deviations  were  as  low  as  2%. 

Bayouth  et.al.^  concluded  in  their  study  that  interleaf  leakage  ranged  from  1 .0%- 
1.5%,  while  intraleaf  transmission  and  leaf-end  transmission  ranged  from  0.6%-0.8%  and 
0.8%-2.7%  respectively. 

Modeling  by  treatment  planning  systems  must  be  able  to  account  for  these 
inherent  small  field  dosimetry  issues.  While  there  has  been  some  work  done  to  describe 
the  large  field  size  dosimetric  characteristics  of  MLC’s^’®’^ ,  there  is  still  much  to  be 
studied  in  order  to  describe  the  small  field  dosimetric  characteristics  of  MLC’s  and 
whether  or  not  significant  differences  exist.  The  objective  of  this  work  was  to  investigate 
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whether  or  not  a  clinically  significant  ditference  exists  in  the  measured  dose  when  using 
the  current  ABAC  Pinnacle^  IMRT  model,  which  is  based  on  large  field  size  data,  versus 
using  a  small  field  models,  based  on  small  field  measurements. 


TT  M  ATERIAT  S  ANtt  METHODS 


A.  EXPOSURES 

Our  institution,  Wilford  Hall  Medical  Center,  Lackland  AFB,  TX  recently 
purchased  a  Varian  (Varian  Medical  Systems,  Palo  Alto,  CA)  21  EX  Linac  with  120  leaf 
MLC IMRT  capability.  Each  leaf  and  carriage  was  independently  controlled.  Leaf 
movement  was  in  the  “X”  direction  (parallel  with  the  lower  jaws).  Maximum  leaf  retract 
position  for  the  Linac  21EX  is  20.1  cm  from  the  beam  centerline  and  maximum  leaf 
extend  position  is  -20.0  cm  over  the  beam  centerline.  The  leaf  thickness  is  5  mm;  the  leaf 
height  is  60.0  mm  with  an  end  radius  of  80.0  mm.  The  leaf  tongue  and  groove  offsets  are 
0.4  mm.  The  coincidence  light  field  vs.  x-ray  is  1  mm.  Penumbra  is  20-80%  leaf  end;  7 
mm  or  less  for  a  10cm  x  10cm  field  for  x-ray  energies  <10MV,  and  8.5  mm  or  less  for 
higher  energies.  The  average  leaf  transmission  is  <2.5%  and  the  maximum  interleaf 
leakage  is  <4.0%.* 

Small  field  dosimetry  data  was  collected  using  Kodak  EDR2  ReadyPak®  film 
(Eastman  Kodak  Company,  Rochester,  NY).  All  exposures  were  made  at  100cm  SAD, 
6MV,  300  monitor  units  (1  MU  =  IcGy),  at  depths  of  1 .5  cm,  5.0  cm,  10  cm,  15cm,  and 
20  cm.  Field  sizes  ranged  from  5  cm^  to  1  cm^  reducing  the  field  size  by  Icm^ 
increments.  Two  sets  of  exposures  were  made,  one  using  primary  jaws  to  collimate,  the 
other  using  the  MLC  to  collimate.  The  films  were  processed  with  a  Kodak  5000  rapid 
film  processor.  A  calibration  film  was  taken  in  conjunction  with  the  dosimetry  films. 
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Film  calibration  was  performed  using  the  method  presented  by  Hanny  et.af. 

Once  processed,  films  were  then  scanned  using  a  VIDAR  DosimetryPRO™ 

Advantage  Film  Digitizer  (Vidar  Systems  Corporation,  Herndon,  VA).  This  scanner  uses 
32-bit  data  path  outputting  to  16  bits  of  greyscale  for  high  resolution.  The  scanned 
images  were  then  captured,  cropped,  and  saved  with  a  75  dots  per  inch  (dpi)  resolution. 
The  net  optical  densities  were  plotted  as  a  function  of  the  measured  dose  to  produce  the 
sensitometric  curve  data.  Data  were  then  analyzed  by  a  nonlinear  least  squares  fit  of  the 
beam  sensitometric  curve  data  to  the  double-double  hit  film  response  model  described  by 
Zhu  et.al.'}'^ 

OD  =  ODi,n,ax  [1-  6-“’°  (1  +  aiD)]  +  OD2,max  [1  “  (1  +  a2D)] 

where  the  OD  is  the  optical  density  of  the  EDR2  film,  ODi^max  and  OD2,max  are  the 
maximum  optical  densities,  aj  and  ai  are  coefficients  related  to  film  sensitivity  (cGy'), 
and  D  is  the  delivered  dose  (cGy). 

With  optical  densities  converted  to  dose,  beam  profiles  were  then  generated.  For 
each  profile  it  was  assumed  that  median  value  of  the  dose  values  equal  to  or  greater  than 
50%  of  the  maximum  dose  represented  the  mid-point  of  the  profile.  This  point  was 
assigned  an  “X”  value  of  zero  and  subsequent  “X”  positions  were  assigned  using  the 
distance  between  each  dot  to  form  the  “X”  axis  of  the  profile.  The  data  was  then 
converted  into  ASCII  format  and  imported  into  ADAC  Pinnacle^. 

Using  this  profile  data,  two  separate  models  were  created: 

1.)  Jaw  model-  Small  field  size  model  constructed  with  data  from  field  sizes  of  5 
cm  X  5  cm  and  smaller  using  primary  jaw  collimation 
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2.)  MLC  model-  Small  field  size  model  constructed  with  data  from  field  sizes  of 
5  cm  X  5  cm  and  smaller  using  MLC  collimation 
A  series  of  comparisons  were  made  using  Pinnacle^®  IMRT  treatment  planning 
software  to  test  the  Jaw  model,  and  MLC  model,  against  the  standard  ABAC  model. 

The  ADAC  model  is  based  on  large  field  size  data,  i.e.,  5  cm  x  5  cm  and  larger.  Even 
though  technically  all  three  models  were  designed  in  ADAC  Pinnacle,  to  avoid  confusion 
they  will  be  labeled  this  point  forward  as  the  Jaw,  MLC,  and  ADAC  models. 

The  following  describes  algorithm  that  ADAC  Pinnacle^®  IMRT  treatment 
planning  software  uses.  The  Pinnacle^®  Collapsed  Cone  Convolution  Superposition  dose 
algorithm  is  based  on  the  work  of  Mackie,  et.al}^  Rather  than  correcting  measured  dose 
distributions,  the  Convolution  Superposition  algorithm  computes  dose  distributions  from 
first  principles  and,  therefore,  can  account  for  the  effects  of  beam  modifiers,  the  external 
patient  contour,  and  tissue  heterogeneities  on  the  dose  distribution.  The  Convolution 
Superposition  dose  model  consists  of  four  parts: 

•  Modeling  the  incident  energy  fluence  as  it  exits  the 
accelerator  head. 

•  Projection  of  this  energy  fluence  through  the  density 
representation  of  a  patient  to  compute  a  TERMA  (Total 
Energy  Released  per  unit  Mass)  volume. 

•  A  three-dimensional  superposition  of  the  TERMA  with  an 
energy  deposition  kernel  using  a  ray-tracing  technique  to 
incorporate  the  effects  of  heterogeneities  on  lateral  scatter. 

•  Electron  contamination  is  modeled  with  an  exponential 
falloff  which  is  added  to  the  dose  distribution  after  the 
photon  dose  is  computed.^ 
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The  series  of  comparisons  looked  at  how  the  Jaw  and  MLC  models  each 
compared  to  the  ADAC  model.  The  first  test  was  to  compare  the  lateral  beam  profiles,  at 
a  depth  of  1 .5  cm  (dmax).  The  Jaw  model  profile  was  compared  to  the  ADAC  model  and 

likewise  the  MLC  model  was  compared  to  the  ADAC  model. 

To  provide  further  analysis  of  the  three  models,  a  single  plane  from  two  separate 
IMRT  plans  (a  head  and  neck  plan  and  a  prostate  plan)  were  used.  For  the  head  and  neck 
IMRT  plan  the  right  lateral  view  was  used  to  generate  three  separate  planar  dose  profiles, 
each  based  on  one  of  the  three  models  (ADAC,  MLC,  Jaw).  For  the  prostate  IMRT  plan 
the  PA  view  was  used  to  do  the  same. 

Using  ADAC  Pinnacle  Smarts  im  to  determine  the  effect  of  the  models  on  the 
conversion  process,  the  calculated  optimized  plans  were  converted  to  calculated 
deliverable  plans  for  each  of  the  models.  The  calculation  parameters  used  in  ADAC 
were;  Error  tolerance  2%,  Minimum  segment  area  2cm^,  Minimum  equivalent  square 
2cm,  Leaf  field  edge  overlap  0.5  cm.  Planar  dose  images  were  then  constructed  to 
determine  if  any  significant  discrepancies  were  visually  appreciated. 

In  order  to  compare  the  calculated  deliverable  plans  for  each  of  the  models,  the 
planar  dose  data  was  imported  into  Dose  Lab  software.  Version  3.05.  The  Jaw  and  MLC 
models’  calculated  deliverable  data  was  compared  to  the  ADAC  model  calculated 
deliverable  data  using  vertical  and  horizontal  dose  profiles,  and  isodose  maps. 

Using  the  planar  dose  data  from  each  of  the  respective  models  and  planes,  MLC 
segment  exposures  were  made  using  the  same  film  and  set  up  as  the  previous  profile 
exposures.  The  films  were  processed,  scanned,  and  sensitometric  data  determined.  To 
determine  how  accurately  each  model  predicted  its’  measured  data,  a  comparison  was 
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made  between  the  calculated  deliverable  data  and  the  deliverable  measured  data.  Using 
Dose  Lab  software  a  NAT  (Normalized  Agreement  Test)  and  DTA  (distance  to 

agreement)  test  were  performed  for  each  model. 

DTA  calculations  are  based  on  the  number  of  contour  lines  assigned.  For  this 
project  a  total  of  7  isodose  lines  were  chosen  (25cGy,  65cGy,  75cGy,  85cGy,  90cGy, 
98cGy).  Coordinates  of  isodose  lines  are  generated  for  a  particular  dose  level  on  both  the 
measured  and  computed  images.  A  point  on  the  measured  isodose  line  is  chosen,  and  the 
distance  fi'om  every  computed  isodose  point  to  the  measured  point  is  calculated.  The 
DTA  image  is  assigned  the  value  of  the  minimum  distance  found  at  the  same  coordinate 
as  the  point  in  the  measured  image.  This  is  repeated  for  all  of  the  measured  isodose 
points  and  other  isodose  lines. 

The  NAT  image  is  based  on  local  percent  difference  and  DTA  values.  If  a  pixel 
passes  the  percent  difference  or  DTA  criteria,  its  NAT  value  is  zero.  To  make  it  more 
useful  clinically,  the  NAT  value  is  assigned  zero  for  cold  areas  outside  the  PTV  (defined 
as  areas  with  computed  doses  less  than  75%  of  the  maximum  dose  as  defined  in  the 
patient  information)  and  the  deviations  from  dose  and  distance  criteria  are  scaled  by 
dose.  The  NAT  index  is  the  average  NAT  value  divided  by  the  average  value  of  the  dose 
scaling  matrix.  This  single  number  represents  the  overall  goodness  of  the  dose  agreement 
(lower  numbers  are  better). 

Also  compared  were  vertical  and  horizontal  dose  profiles,  as  well  as  isodose 
maps.  For  completeness  a  2D  histogram  of  the  percent  of  pixels  that  had  a  3%  dose  or  a 
3mm  distance  difference  was  generated.  While  it  was  helpful  in  determining  dose  and 
distance  criteria,  its  loss  of  spatial  information  should  be  considered. 


in.  RESULTS 


A.  JAW  PROFILES 

The  following  figures  (1-5)  are  lateral  beam  profiles  generated  by  the  measured 
data  using  jaws  for  collimation,  and  compared  to  the  AD  AC  model.  Figures  6  and  7 
represent  depth  dose  data  comparisons  for  the  2  cm^  and  5  cm^  field  sizes. 

For  the  smaller  field  sizes,  1  cm^  -  3  cm^  the  measured  data  from  jaw  collimation 
varied  from  AD  AC  by  a  range  of  2%  -  1 1%.  For  the  I  cm^  profile  he  measured  data 
reflected  a  higher  dose  in  the  central  part  of  the  field  but  a  lower  dose  in  the  penumbra 
region.  For  the  2  cm^  field  size,  there  was  good  agreement  in  the  central  part  of  the 
profile  with  slight  deviation  at  the  edges.  The  3  cm^  profile  reflected  the  greatest 
inconsistency  the  measured  data  reflecting  a  lower  dose  profile  compared  to  the  AD  AC 
model.  For  the  4  cm^  and  5  cm^  field  sizes,  the  measured  profiles  compared  well  with  the 
AD  AC  model  with  differences  =/<  2%.  It  was  expected  that  for  these  larger  field  sizes, 
there  would  be  good  agreement.  The  differences  seen  with  the  smaller  field  sizes  suggest 
that  AD  AC  predicts  a  higher  dose  than  was  actually  measured  with  film  at  the  central 
50%  of  the  beam.  The  depth  dose  comparison  shows  good  agreement  between  the  Jaw 
and  AD  AC  model. 
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B.  MLC  PROFILES 

The  following  figures  (8-12)  are  beam  profiles  generated  by  the  measured 
data  using  the  MLC  for  collimation,  and  compared  to  the  ABAC  model.  Figures  13  and 
14  represent  depth  dose  data  comparisons  for  the  2  cm^  and  5  cm^  field  sizes. 

For  the  1  cm^  field  size  the  measured  data  ifom  MLC  collimation  varied  from 
ABAC  by  a  range  of  2%  -  12%  along  the  central  50%  of  the  beam.  The  MLC  dose 
profile  was  consistently  greater  above  the  50%  maximum  dose  point.  Below  that  dose 
point  ABAC  predicted  a  narrower  penumbra  than  actually  measured.  For  2  cm  and  3 
cm^  field  sizes  there  was  good  agreement  between  ABAC  and  MLC  data  with  a 
difference  of  2%-4%  overall.  For  the  4  cm^  and  5  cm^  it  was  interesting  to  see  that  the 
MLC  data  reflected  a  profile  consistently  less  than  the  ABAC  model  predicted.  The 
depth  dose  comparison  shows  good  agreement  between  the  MLC  and  ABAC  model. 
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C.  HEAD  AND  NECK  PLANAR  DOSE  DATA  AND  IMAGES 
1.  COMPARISON  OF  PLANAR  IMAGES 
The  following  figures  (15-17)  are  planar  images  of  a  right  lateral  neck  view  from 
an  actual  head  and  neck  IMRT  plan.  The  first  image  represents  the  AD  AC  model.  The 
following  images  reflect  the  same  view  but  using  the  Jaw  model  and  MLC  model 
respectively. 
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2.  COMPARISON  OF  CALCULATED  DELIVERABLE  DATA 
Figure  18  represents  a  vertical  dose  profile  comparison  of  the  calculated 
deliverable  data  using  the  ABAC  model  versus  the  Jaw  model. 
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3.  COMPARISON  OF  CALCULATED  DELIVERABLE  VS. 

DELIVERABLE  MEASURED  DATA 
To  determine  which  model  best  predicted  its  measured  dose  from  film,  the 
calculated  deliverable  data  was  compared  to  the  deliverable  measured  data  for  each 
model.  The  following  figures  represent  this  series  of  comparisons.  All  compare  the  same 
right  lateral  neck  IMRT  view. 

Figures  24-29  represent  the  tests  used  to  compare  the  AD  AC  model 
calculated  deliverable  plan  versus  the  ABAC  model  deliverable  measured  data  from  film. 
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When  comparing  the  AD  AC  calculated  deliverable  plan  versus  the  ADAC 
deliverable  measured  data,  the  DTA  scores  were  only  in  the  50%  range.  There  is 
noticeable  dose  difference  noted  in  the  left  side  of  the  DTA  and  NAT  images.  When 
viewing  the  2D  dose  histogram,  there  are  a  significant  amount  of  pixels  which  fall  outside 
the  thresholds  of  3%  or  3mm. 

Figures  30-35  represent  the  comparison  between  the  MLC  calculated 
deliverable  plan  versus  the  MLC  deliverable  measured  data  from  film. 
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With  the  distance  to  agreement  scores  in  the  70%  range,  the  MLC  model  was  a 
better  predictor  for  its  measured  data  than  the  AD  AC  model  was  for  its  measured  data. 

Figures  36-41  represent  the  comparison  between  the  Jaw  model  calculated 


deliverable  dose  versus  the  deliverable  measured  data. 
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The  distance  to  agreement  scores  were  as  high  as  the  MLC  comparison.  The  Jaw 
model  was  a  better  predictor  for  its  measured  data  than  the  AD  AC  model  was  for  its 
measured  data.  The  following  summarizes  the  NAT  and  DTA  scores  for  the  models. 


Calculated  deliverable  vs.  deliverable  measured  from  film 


D.  PROSTATE  PLANAR  DOSE  DATA  AND  IMAGES 
1.  COMPARISON  OF  PLANAR  IMAGES 
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The  following  figures  (43-45)  are  planar  images  of  a  posterior  to  anterior 
view  from  an  actual  prostate  IMRT  plan.  The  first  image  represents  the  AD  AC 
model.  The  following  images  reflect  the  same  view  but  using  the  Jaw  model  and 
MLC  model  respectively. 
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The  images  visually  reflect  good  agreement  between  all  three  models  used. 
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2.  COMPARISON  OF  CALCULATED  DELIVERABLE  DATA 
Figure  46  represents  a  vertical  dose  profile  comparison  of  the  calculated 
deliverable  data  using  the  AD  AC  model  versus  the  Jaw  model. 
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3.  COMPARISON  OF  CALCULATED  DELIVERABLE  VS. 

DELIVERABLE  MEASURED  DATA 
To  determine  which  model  best  predicted  its’  measured  dose  the  calculated 
deliverable  data  was  compared  to  the  deliverable  measured  data  for  each  model.  The 
following  figures  represent  this  series  of  comparisons.  All  compare  the  same  prostate 
posterior  to  anterior  IMRT  view. 

Figures  52-57  represent  various  tests  used  to  compare  the  ADAC 
calculated  deliverable  plan  versus  the  deliverable  measured  data  from  the  ADAC  model. 
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For  the  prostate  plan,  the  AD  AC  optimized  plan  did  much  better  at  predicting  the 
measured  dose  versus  the  neck  plan.  The  distance  to  agreement  score  was  79%.  Looking 
at  the  2D  histogram,  there  were  some  borderline  dose  discrepancies,  which  fell  outside 
the  distance  threshold  of  3mm. 

Figures  58-63  represent  the  comparison  between  the  MLC  model  calculated 
deliverable  dose  versus  the  deliverable  measured  dose. 
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The  high  distance  to  agreement  score  of  97%  for  this  comparison  suggests  that  the 
MLC  model  was  a  very  accurate  predictor  for  measured  dose.  This  is  significantly  higher 
than  the  AD  AC  model  scored. 


Figures  64-69  represent  the  comparison  between  the  Jaw  planar  dose  data  and  the 


Jaw  measured  dose  data. 
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The  Jaw  model  had  a  distance  to  agreement  score  slightly  less  than  the  MLC 
model  but  still  predicted  measured  dose  from  its  own  plan  very  well.  This  was  also 
significantly  higher  than  the  ADAC  model  scored.  The  following  summarizes  the  NAT 
and  DTA  scores. 


Calculated  deliverable  vs.  deliverable  measured  from  film 


TV.  DISCUSSION 


Historically,  treatment  planning  systems  do  not  include  beam  data  for  field  sizes 
less  than  3cm  x  3cm  in  the  beam  model  because  it  is  rare  to  plan  with  such  small  fields. 
However,  in  IMRT,  major  portions  of  the  treatment  can  be  given  with  fields  as  small  as 
1cm  X  ]  cm.  These  small  field  sizes  and  their  inherent  dosimetric  peculiarities  are 
relatively  new  to  treatment  planning  algorithms  and  have  to  be  carefully  addressed  in  the 
beam  model  in  order  to  get  accurate  dose  calculations  for  IMRT. 

The  method  described  in  this  report  was  a  method  used  at  our  institution  to  test 
the  ADAC  Pinnacle^  Model  (AD AC),  a  model  which  we  implemented  based  on  large 
field  dosimetry,  against  two  other  models  (MLC  and  Jaw),  based  on  small  field  dosimetry. 

When  comparing  the  beam  profiles,  there  were  differences  noted  between  the 
three  models  that  ranged  from  2%  -  20%.  This  was  especially  prominent  for  the  smaller 
field  sizes,  i.e.,  Icm^  -  3cm^. 

Using  a  Normalized  Agreement  Test  and  a  Distance  to  Agreement  Test,  we  were 
able  to  determine  that  the  small  field  models  were  better  predictors  for  delivered  dose. 
Although  it  is  ultimately  up  to  the  physicist  and  oncologist  to  deteraiine  what  dose  and 
distance  thresholds  are  acceptable,  by  convention  a  3%  dose  or  3mm  variance  are 
considered  the  norm.’^’’"^  The  Normalized  Agreement  Test  and  Distance  to  Agreement 
scores  were,  on  average,  24%  higher  for  the  Jaw  and  MLC  models  versus  the  ADAC 
model.  These  results  were  expected  since  the  ADAC  model  does  not  include  dosimetry 


67 


68 


data  for  small  field  sizes.  The  fact  that  the  Jaw  model  scored  almost  as  well  as  the  MLC 
model  suggests  that  it  is  the  small  field  modeling  itself  that  contributes  to  the  overall  dose 
and  distance  agreement  and  that  an  MLC  model  versus  a  Jaw  model  provides  only  a 
slightly  better  performance. 

All  three  models  scored  lower  for  the  head  and  neck  plan  versus  the  prostate  plan. 
This  is  likely  a  function  of  the  increased  modulation  for  this  particular  plan.  The  fact  that 
there  was  a  greater  difference  in  scores  with  the  head  and  neck  plan  versus  the  prostate 
plan,  (30%  difference  for  the  head  and  neck  plan  and  17%  difference  for  the  prostate 
plan),  suggests  that  with  a  plan  that  has  increased  modulation,  the  beam  modeling 
becomes  even  more  critical. 

Based  on  our  results,  we  suggest  that  treatment  planning  system  users  evaluate  the 
performance  of  their  beam  models  and  consider  whether  to  incorporate  small  field  data 
into  the  model.  This  is  especially  warranted  with  IMRT  plans  that  have  extensive 
modulation.  Due  to  the  increased  numbers  of  small  field  sizes  generated  by  IMRT  plans, 
the  treatment  planning  system  will  need  small  field  dosimeteric  data  to  accurately  predict 


dose. 


TV.  CONCLUSIONS 


The  objective  of  this  work  was  to  test  whether  ADAC  Pinnacle^  IMRT 
treatment  planning  software  can  accurately  predict  doses  for  small  field  sizes.  It  was  our 
hypothesis  that  the  ADAC  model  would  not  accurately  predict  dose  for  small  fields.  It 
was  fixrther  hypothesized  that  this  inaccuracy  would  be  clinically  significant. 

Our  results  show  that  small  field  models,  using  small  field  dosimetry  data  fi'om  either  jaw 
or  MLC  collimation,  were  better  predictors  for  IMRT  dose  than  the  ADAC  Pinnacle 
model. 

The  data  suggest  that  a  clinically  significant  difference  can  arise  between  IMRT 
plans  that  are  based  on  large  and  small  field  models.  Further  studies  would  be  required  to 
determine  whether  or  not  patient  outcome  would  be  affected  by  the  choice  of  the  beam 
sizes  used  to  model  the  radiation  characteristics. 
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APPENDIX 


The  ABAC  Pinnacle3  Collapsed  Cone  Convolution 
Superposition  Dose  Model 

Todd  McNutt,  Ph.D.  -  Director  of  Product  Development 

ADAC’s  Pinnacles  3D  treatment  planning  system  uses  a  Collapsed  Cone  Convolution 
Superposition  computation  to  determine  the  dose  distribution  in  patients  from  external 
photon  beams.  The  PinnacleS  Convolution  Superposition  dose  model  is  a  true  three 
dimensional  dose  computation  which  intrinsically  handles  the  effects  of  patient 
heterogeneities  on  both  primary  and  secondary  scattered  radiation.  This  computation 
method  is  uniquely  able  to  account  for  dose  distributions  in  areas  where  the  electronic 
equilibrium  is  perturbed,  such  as  tissue-air  interfaces  and  tissue-bone  interfaces.  While 
other  convolution  techniques  account  for  the  effects  of  patient  heterogeneities  on  primary 
radiation,  they  neglect  the  effects  of  heterogeneities  on  scattered  radiation  in  the  final 
dose  distribution.  In  addition,  the  nature  of  the  Convolution  Superposition  dose  model 
makes  it  ideally  suited  to  handle  optimization  and  intensity  modulated  radiation  therapy 
planning. 

The  Convolution  Superposition  Dose  Model 

The  Pinnacles  Convolution  Superposition  dose  algorithm  is  based  on  the  work  of 
Mackie,  et  al.  Rather  than  correcting  measured  dose  distributions,  the  Convolution 
Superposition  algorithm  computes  dose  distributions  from  first  principles  and,  therefore, 
can  account  for  the  effects  of  beam  modifiers,  the  external  patient  contour,  and  tissue 
heterogeneities  on  the  dose  distribution. 
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The  Convolution  Superposition  dose  model  consists  of  four  parts; 

•  Modeling  the  incident  energy  fluence  as  it  exits  the  accelerator  head. 

•  Projection  of  this  energy  fluence  through  the  density  representation  of  a  patient  to 
compute  a  TERMA  (Total  Energy  Released  per  unit  Mass)  volume. 

•  A  three-dimensional  superposition  of  the  TERMA  with  an  energy  deposition  kernel 
using  a  ray-tracing  technique  to  incorporate  the  effects  of  heterogeneities  on  lateral 
scatter. 

•  Electron  contamination  is  modeled  with  an  exponential  falloff  which  is  added  to  the 
dose  distribution  after  the  photon  dose  is  computed. 

The  following  sections  describe  each  part  of  the  model  in  more  detail. 

Modeling  the  Incident  Energy  Fluence  as  it  exits  the  Accelerator 
The  incident  energy  fluence  distribution  is  modeled  as  a  two-dimensional  array  which 
describes  the  radiation  exiting  the  head  of  the  linear  accelerator.  The  parameters  defining 
this  array  is  defined  during  physics  data  modeling.  The  starting  point  for  photon 
modeling  is  a  uniform  plane  of  energy  fluence  describing  the  intensity  of  the  radiation 
exiting  the  accelerator  head.  The  fluence  model  is  then  adjusted  to  account  for  the 
flattening  filter,  the  accelerator  head,  and  beam  modifiers  such  as  blocks,  wedges  and 
compensators. 

•  The  “horns”  in  the  beam  produced  by  the  flattening  filter  are  modeled  by  removing  an 
inverted  cone  from  the  distribution. 

•  Off-focus  scatter  produced  in  the  accelerator  head  is  modeled  by  defining  a  2D 
Gaussian  function  as  a  scatter  source  and  adjusting  the  incident  energy  fluence  based  on 
the  portion  of  the  Gaussian  distribution  visible  from  each  point  in  the  incident  energy 
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fluence  plane. 

•  The  geometric  penumbra  is  modeled  by  convolving  the  array  with  a  focal  spot 
blurring  function. 

•  During  planning,  the  shape  of  the  field  produced  by  blocks  or  multi-leaf  collimators 
is  cut  out  of  the  array  leaving  behind  the  corresponding  transmission  through  the 
shape-defining  entity. 

•  Beam  modifiers  such  as  wedges  and  compensators  are  included  in  the  array  by 
attenuating  the  energy  fluence  by  the  corresponding  thickness  of  the  modifier.  For  static 
wedges  and  compensators,  a  radiological  depth  array  is  also  stored  which  allows  for 
proper  modeling  of  the  beam  hardening  due  to  the  presence  of  the  beam  modifiers  during 
the  projection  of  the  incident  fluence  array. 

Dynamic  beam  delivery  with  intensity-modulation  or  dynamic  wedges  is  easily  handled 
using  the  incident  energy  fluence  array.  For  these  beams,  the  radiological  depth  array  is 
not  needed  to  account  for  beam  hardening. 

Projection  of  Energy  Fluence  through  a  CT  Patient  Representation 
The  incident  energy  fluence  plane  is  projected  tlirough  the  CT  patient  representation  and 
attenuated  using  mass  attenuation  coefficients.  These  coefficients  are  stored  in  a 
threedimensional  lookup  table  as  a  function  of  density,  radiological  depth,  and  off-axis 
angle.  Patient  heterogeneities  are  taken  into  account  with  the  density  dependence.  Beam 
hardening  through  the  patient  is  accounted  for  with  the  radiological  depth  dependence, 
and  the  off-axis  softening  of  the  energy  spectrum  is  produced  with  the  off-axis  angle 
dependence.  To  account  for  the  changes  in  the  photon  energy  spectrum  at  different 
locations  in  the  beam,  the  mass  attenuation  coefficient  lookup  table  is  produced  using 


72 


a  weighted  sum  of  several  mono-energetic  tables.  The  TERMA  (Total  Energy  Released 
per  unit  Mass)  volume  is  computed  by  projecting  the  incident  energy  fluence  through  the 
patient  density  volume  using  a  ray-tracing  technique.  A  given  ray’s  direction  is 
determined  based  on  the  position  of  the  radiation  source  and  the  particular  location  in  the 
incident  fluence  plane.  At  each  voxel  in  the  ray  path,  the  TERMA  is  computed  using  the 
attenuated  energy  fluence  along  the  ray  and  the  mass  attenuation  coefficient  at  the 
particular  density,  radiological  depth,  and  off-axis  angle. 

3D  Superposition  of  an  Energy  Deposition  Kernel 

The  three-dimensional  dose  distribution  in  the  patient  is  computed  by  superposition  of  the 
TERMA  volume  with  the  energy  deposition  kernel.  The  kernel  represents  the  spread  of 
energy  from  the  primary  photon  interaction  site  throughout  the  associated  volume.  Poly- 
energetic  kernels  are  produced  by  combining  a  series  of  Monte  Carlo-generated  mono- 
energetic  energy  deposition  kernels.  The  superposition  is  carried  out  using  a  ray  tracing 
technique  similar  to  that  used  in  the  projection  of  the  incident  enei'gy  fluence.  The  kernel 
is  inverted  so  that  the  dose  can  be  computed  in  only  a  portion  of  the  patient  (TERMA) 
volume  if  desired.  This  allows  for  point  dose  computation  and  decreases  computation 
time.  The  rays  from  the  dose  deposition  site  are  cast  in  three  dimensions.  At  each  voxel 
of  the  TERMA  traversed  along  a  ray,  the  contribution  of  dose  to  the  dose  deposition  site 
is  computed  and  accumulated  using  the  TERMA  and  the  kernel  value  at  the  current 
radiological  distance.  Using  the  radiological  distance  along  the  ray  also  allows  the  kernel 
to  be  scaled  to  account  for  the  presence  of  heterogeneities  with  respect  to 
scattered  radiation  in  all  directions.  The  dose  computation  described  above  determines 
the  dose  from  a  single  beam.  Multiple  beams  are  computed  independently  and  the  entire 
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3D  dose  distribution  is  created  by  adding  the  dose  from  each  beam  together  according  to 
the  corresponding  beam  weight. 

Adaptive  Convolution  Superposition 

An  Adaptive  Convolution  Superposition  approach  has  also  been  implemented  in 
Pinnacles.  This  uses  the  calculation  technique  described  above  with  some  slight 
modifications.  The  speed  of  the  computation  is  increased  by  adaptively  varying  the 
resolution  of  the  dose  computation  grid  depending  on  the  curvature  of  the  TERMA  and 
dose  distribution.  First,  the  dose  in  a  coarse  3D  grid  is  computed  and  then  the  curvature  in 
the  TERMA  distribution  is  assessed.  In  regions  where  the  curvature  is  high,  the  dose  is 
computed  at  intermediate  points  to  provide  higher  resolution.  The  system  adaptively 
increases  the  resolution  in  regions  of  high  curvature  until  an  acceptable  resolution  is  used. 
In  regions  of  low  curvature,  the  dose  is  interpolated  from  the  coarse  dose  grid.  This 
technique  decreases  the  computation  time  by  a  factor  of  2-3  without  compromising  the 
accuracy  of  the  Convolution  Superposition  calculation  in  the  presence  of  heterogeneities. 

Other  Model-Based  Algorithms 

Other  model-based  algorithms,  including  3D  Fast  Fourier  Transform  (FFT)  techniques 
and  differential  pencil  beam  models  which  use  FFTs  on  two-dimensional  planes,  use  a 
projection  of  the  incident  energy  fluence  similar  to  that  used  in  Pinnacle3  to  determine 
the  TERMA  volume.  They  differ  in  that  they  do  not  use  the  superposition  technique  in 
the  convolution  process.  The  FFT  techniques  require  the  assumption  of  an  invariant 
kernel,  which  inherently  assumes  a  homogeneous  density  representation  during  the 
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convolution  process.  This  technique  reduces  the  accuracy  of  the  computation  because  it 
ignores  the  effects  of  heterogeneities  on  laterally  scattered  radiation.  Some  post 
computation  corrections  may  be  performed  to  help  alleviate  the  error.  In  contrast,  by 
scaling  the  rays  from  the  primary  dose  deposition  site,  the  Convolution  Superposition 
method  accurately  and  intrinsically  models  the  effects  of  lateral  scatter  from  tissue 
heterogeneities,  a  requirement  for  calculating  dose  from  conformal  and  intensity 
modulated  fields. 

AD  AC  Laboratories 
540  Alder  Dr. 

Milpitas,  CA  95035 
Tel:  (408)321-9100  (800)538-8531 
Fax:  (408)577-0907  www.adaclabs.com 
MBA-XOOlO,  Rev.  B 

Although  the  FFT  algorithms  have  a  fast  computation  speed  per  computation  point,  they 
require  full  computation  of  dose  over  the  entire  TERMA  volume.  For  irregular  field 
calculations,  point  doses,  plane  doses,  or  other  situations  where  the  planner  is  only 
interested  in  the  dose  to  smaller  regions,  the  FFT  algorithms  still  require  the  dose  to  be 
computed  over  the  entire  volume.  The  Convolution  Superposition  method  can  accurately 
compute  the  dose  to  a  single  point  with  the  inverted  energy  deposition  kernel.  Therefore, 
depending  on  the  planning  situation,  the  desired  dose  calculation  may  be  faster  using  the 
Convolution  Superposition  calculation  than  when  using  the  FFT  calculation. 

The  ability  to  define  a  smaller  calculation  matrix  also  results  in  a  calculation  speed 
advantage  with  the  Convolution  Superposition  model  for  inverse  planning  and 
optimization  of  intensity  modulated  beams  where  computation  of  dose  need  only 
be  performed  in  limited  regions  during  the  optimization  process. 
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